Fendley, Ludwig and Saleur reply to Skorik's com- 
ment: 

In a series of papers Q , we have used a novel approach, 
combining the Bethe ansatz with a kinetic (Boltzmann) 
equation, in order to compute exactly transport proper- 
ties in a Luttinger liquid with an impurity (also known as 
the boundary sine-Gordon model [BSGM]). In a recent 
comment Skorik claimed that while our results were 
correct in linear response, there was a "serious flaw" in 
our calculation at finite voltage. We explain here why 
Skorik's arguments are inappropriate, and that there is 
no flaw in our work. 

In our approach, the effect of the voltage corresponds 
to injecting into the sample left movers and right movers 
from different reservoirs, at chemical potentials, hl,r — 
respectively! |,[| . We use a Landauer-Biittiker scat- 
tering approach || . By a non-linear change of basis called 
"folding" §, we map this to a problem with reservoirs 
setting different chemical potentials for the quasiparticles 
of charge ±. These quasiparticles are || the "massless 
limit" of the usual sine-Gordon soliton and antisoliton 
states. Because the model is integrable, these scatter 
one-by-one off each other and off the impurity, conserv- 
ing their momentum and either conserving or switching 
charge, with exactly known reflection and transmission 
S'-matrices. This permits the construction of the exact 
(scattering) eigenstates of the interacting Hamiltonian, 
describing the Luttinger liquid leads plus the impurity, 
in a way analogous to ordinary potential scattering. To 
compute the conductance, we use a Boltzmann equation, 
counting how much charge is transported through the 
impurity in the presence of different populations of ± 
quasiparticles, set by the bias. This is natural, in spite 
of the interacting nature of the problem, because of the 
very simple nature of collisions that follows from integra- 
bility. In fact, we note that subsequent to our original 
work, the Boltzmann equation, in the case of linear re- 
sponse, was derived Q directly from the Kubo formula 
combined with form-factors, i.e. matrix elements of cur- 
rent operators in the quasi particle basis, and complete 
agreement with the earlier results jl| was found. 

The use of the Boltzmann equation per se was not 
criticized by Skorik. As far as we understand it, his con- 
cern is that in the Boltzmann equation we used scatter- 
ing matrix elements "at zero voltage" , whereas, for finite 
voltage, the S matrix itself might potentially acquire a 
voltage dependence, determined by the V dependent fill- 
ing of the ground state. We will show that this is in 
fact not the case. First, we note that following Skorik's 
logic, the same criticism could seemingly be made for 
the linear-response calculation at non-zero temperature. 
Here, the ground state is not shifted by the voltage, but 
physical properties do not depend so much on the zero- 
temperature ground state as on the states whose filling 
fractions are the ones of thermal equilibrium (sec Q for 
more details on this). Nevertheless, Skorik agrees that 
the zero-temperature 5* matrix is the appropriate one to 



use here, and that it does not acquire any sort of tem- 
perature dependence. As discussed in detail previously 
jl| , the only effect of the temperature on the Boltzmann 
equation is the appearance of the non-trivial filling frac- 
tions. To use another language, we are still describing 
the problem in terms of the "bare" particles (we mean 
bare in the sense of thermal fluctuations, not quantum 
fluctuations). The integrability means that these bare 
particles are a valid basis for the problem and that they 
still scatter one-by-one off each other and off the impu- 
rity despite the macroscopic number of particles excited 
around a given particle at non-zero temperature. (More 
precisely, this means that we do calculations in the ex- 
tremely dilute limit where a particle description is valid; 
the central assumption of this and all thermodynamic 
Bethe ansatz computations is that no phase transition in- 
terferes with the continuation of the result to the regime 
of finite densities.) 

Consider now a non-zero voltage, and let us discuss the 
simplest case of T = 0. The Fermi sea is, indeed, filled 
with bare (with respect to the zero-voltage sea) quasi- 
particles, as discussed by Skorik. However, each of these 
scatters with the same, field-independent S'-matrix on 
the impurity. It is this S matrix that is used to build the 
asymptotic states, it is therefore the one that appears in 
the Boltzmann equation in 0. Of course, one can be 
interested in the S-matrix for scattering particles excited 
on top of the finite-field sea: by scattering these particles 
through the impurity and the sea, one gets then a field 
dependent result. A computation of that type was per- 
formed for the the Kondo problem in a magnetic field in 
|Iof , as mentioned by Skorik. Such a dressed S matrix has 
actually appeared in other computations we made, such 
as low frequency AC properties [||. But as far as the 
DC conductance is concerned, the only object necessary 
is the one to build asymptotic states, and it is the "bare 
S matrix" as used in jj]. In fact, even if one wanted to 
use the "dressed" S matrix in our Boltzmann equation, 
nothing would change. This is because, at the reflection 
less points of the sine- Gordon model, to which we have 
restricted, the bulk scattering is diagonal, so the dress- 
ing is a mere phase (it is also the phase that is discussed 
in fLOfl .) Since in our Boltzmann-equation approach the 
DC current depends only on probabilities, i.e. modulus 
squared of S matrix elements, the Independence of this 
S'-matrix would drop out anyway. This, we think, inval- 
idates Skorik's criticism. 

There is one major difference between non-zero tem- 
perature and non-zero voltage, which might have created 
some confusion: As emphasized correctly by Skorik, di- 
agonalizing the BSGM directly at nonzero voltage is a 
difficult exercise, since the boundary interaction does not 
conserve the charge. It was not done in jij: there, we ar- 
gued instead that, physically, the role of the voltage is 
to fix the populations of quasiparticles; it is applied far 
from the impurity, and hence, since the excitations are 
localized solitons, in a region where the effect of the im- 
purity is negligible. Thus at least when computing DC 



transport properties through the impurity, the effect of 
the voltage is to shift the Fermi sea but does not interfere 
with the one-by-one scattering off the impurity (this is 
actually one of the basic ideas underlying the Landauer- 
Biittiker approach). 

In fact, it is also possible to diagonalize the BSGM 
with a voltage, and to confirm this physical argument 
|^l|,|U^Jl^] . The main idea is to observe that the phys- 
ical properties of the BSGM are the same as the prop- 
erties of another model, where the boundary interaction 
cos 4> is replaced by another interaction of the Kondo type 
S + e~ l ^ + S~e 1 ^, where the spin is taken in, either a cyclic 
representation of the U q sl(2) algebra Jljj, or a represen- 
tation of the oscillator algebra jL2j. The key property 
of these representations is that all monomials (of total 
vanishing charge) in S have the same trace. One can 
then compare the perturbative computations of proper- 
ties like the partition function, or the conductance us- 
ing the Keldysh formalism juj in the two models, and 
prove their equivalence (this method is used in jl5) in 
the case g = i , where the are equivalent to boundary 
fermions). The advantage of this new formulation is that 
there is a conserved charge, the sum of the quasiparticles 
charge and the spin S z of the boundary degree of free- 
dom. The voltage is then included by shifting the field 
4> — > <f> + QVt in the boundary interaction |p Ul4|Jiq| . As 
discussed in the appendix of M , this is equivalent to not 
shifting </>, but applying a magnetic field on the boundary 
spin, h = gV . Still in the appendix of 0, it is explained 
how the problem with a field applied only to the bound- 
ary spin, and the problem with a field applied to the 
spin and another field applied to the bulk U(\) charge 
are related by a unitary transformation (this is closely 
related to the behavior of the electrons-impurity suscep- 
tibility in the Kondo problem, as discussed recently in 
[ p~6|] . This transformation shifts the overall charge by a 
constant, but does not modify difference of C7(l) charges). 
We can then consider the problem where the field is cou- 
pled to the conserved charge, which is of course trivial 
to diagonalize. Asymptotic states can then be written 
explicitly for this auxiliary problem with impurity and 
voltage; they involve of course the same S matrix as the 
ones without voltage (the same occurs in the g — \ case 
p5|). Scattering eigenstates in the traditional sense can 
be constructed, and the Landauer-Biittiker approach can 
then be applied to compute the current, the DC fluctua- 
tions |l7j] , and some of the AC properties i_B . 

In conclusion, the formula proposed in does not suf- 
fer from Skorik's criticisms. For the reader who does not 
want to follow the detailed arguments presented above, 
we observe that our formula has passed successfully more 
tests than recognized in Skorik's comment. In addition to 
the case free-fermion case g = h, our formula reproduces 
the correct result for g — ► and g = 1. Even though the 
latter case consists of free fermions in the unfolded ver- 
sion of the problem, it is highly non trivial in our folded 
point of view: indeed, the solitons and antisolitons scat- 
ter then with an SU(2) invariant S matrix, identical to 



the one in the Kondo problem. The former case, while 
simple to study in the classical limit, is highly non-trivial 
from the point of view of integrability: it involves taking 
the complicated limit of an interacting quantum problem, 
where an infinity of quasiparticles scatter, all with non 
trivial S matrices, and the filling fractions of the solitons 
and antisolitons are far from simple. The final result for 
the g — ► limit, expressed in terms of Bessel functions of 
imaginary arguments) ]l^ , p^[ , agrees with the result ob- 
tained using a Fokker-Planck equation jl9| . In addition, 
for any g our formula displays the right behavior in the 
strong and weak backscattering limits, the existence of 
a maximum for large enough voltage which is expected 
on physical grounds H. Finally, it also obeys a duality 
(proven at T = in and related with very natural 
analyticity conjectures at T > in |12|) between weak 
and strong backscattering, where the appropriate tun- 
neling particles are Laughlin quasiparticles and electrons 
respectively. This duality, while not established rigor- 
ously prior to our work, is considered highly desirable on 
physical grounds 14 lq]. We thus see no reason to cast 
doubts on our result, and hope that sooner or later, it 
will be compared favorably with numerical simulation or 
experimental data. 
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